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Abstract 



In these lectures I will give an introduction to Feynman integrals. In the first part of the 
course I review the basics of the perturbative expansion in quantum field theories. In the 
ly-^ . second part of the course I will discuss more advanced topics: Mathematical aspects of loop 

O I integrals related to periods, shuffle algebras and multiple polylogarithms are covered as well 

as practical algorithms for evaluating Feynman integrals. 



1 Introduction 

In these lectures I will give an introduction to perturbation theory and Feynman integrals occur- 
ring in quantum field theory. But before embarking onto a journey of integration and special 
function theory, it is worth recalling the motivation for such an effort. 

High-energy physics is successfully described by the Standard Model. The term "Standard 
Model" has become a synonym for a quantum field theory based on the gauge group SU{3) (8> 
SU{2) i^U{l). At high energies all coupling constants are small and perturbation theory is a 
valuable tool to obtain predictions from the theory. For the Standard Model there are three 
coupling constants, ^i, g2 and ^3, corresponding to the gauge groups U{1), SU{2) and SU{3), 
respectively. As all methods which will be discussed below, do not depend on the specific nature 
of these gauge groups and are even applicable to extensions of the Standard Model (like super- 
symmetry), I will just talk about a single expansion in a single coupling constant. All observable 
quantities are taken as a power series expansion in the coupling constant, and calculated order 
by order in perturbation theory. 

Over the years particle physics has become a field where precision measurements have be- 
come possible. Of course, the increase in experimental precision has to be matched with more 
accurate calculations from the theoretical side. This is the "raison d'etre" for loop calculations: 
A higher accuracy is reached by including more terms in the perturbative expansion. There is 
even an additional "bonus" we get from loop calculations: Inside the loops we have to take into 
account all particles which could possibly circle there, even the ones which are too heavy to be 
produced directly in an experiment. Therefore loop calculations in combination with precision 
measurements allow us to extend the range of sensitivity of experiments from the region which is 
directly accessible towards the range of heavier particles which manifest themselves only through 
quantum corrections. As an example, the mass of top quark has been predicted before the dis- 
covery of the top quark from the loop corrections to electro- weak precision experiments. The 
same experiments predict currently a range for the mass of the yet undiscovered Higgs boson. 

It is generally believed that a perturbative series is only an asymptotic series, which will 
diverge, if more and more terms beyond a certain order are included. However this shall be of no 
concern to us here. We content ourselves to the first few terms in the perturbative expansion with 
the implicit assumption, that the point where the power series starts to diverge is far beyond our 
computational abilities. In fact, our computational abilities are rather limited. The complexity of 
a calculation increases obviously with the number of loops, but also with the number of external 
particles or the number of non-zero internal masses associated to propagators. To give an idea 
of the state of the art, specific quantities which are just pure numbers have been computed up to 
an impressive fourth or third order. Examples are the calculation of the 4-loop contribution to 
the QCD (3-function [1], the calculation of the anomalous magnetic moment of the electron up to 
three loops [2], and the calculation of the ratio 

o(e+e^ — )> hadrons) 

of the total cross section for hadron production to the total cross section for the production of 
a 1-1^1-1 pair in electron-positron annihilation to order O (^3) (also involving a three loop cal- 



culation) [3]. Quantities which depend on a single variable are known at the best to the third 
order. Outstanding examples are the computation of the three-loop Altarelli-Parisi splitting func- 
tions [4, 5] or the calculation of the two-loop amplitudes for the most interesting 2 — )■ 2 pro- 
cesses [6-16]. The complexity of a two-loop computation increases, if the result depends on 
more than one variable. An example for a two-loop calculation whose result depends on two 
variables is the computation of the two-loop amplitudes for e+e^ — )> 3 jets [17-27]. 

On the other hand is the mathematics encountered in these calculations of interest in its own 
right and has led in the last years to a fruitful interplay between mathematicians and physicists. 
Examples are the relation of Feynman integrals to periods, mixed Hodge structures and motives, 
as well as the occurrence of certain transcendental constants in the result of a calculation [28^7] . 
Typical transcendental constants which occur in the final results are multiple zeta values. They 
are obtained from multiple poly logarithms at special values of the arguments. I will discuss these 
functions in detail in these lectures. 

The outline of this course is as follows: In section 2 I review the basics of perturbative quan- 
tum field theory and I give a brief outline how Feynman rules are derived from the Lagrangian 
of the theory. Issues related to the regularisation of otherwise divergent integrals are treated in 
section 3. Section 4 is devoted to basic techniques, which allow us to exchange the integrals 
over the loop momenta against integrals over Feynman parameters. In sect. 5 I discuss how the 
Feynman parametrisation for a generic scalar /-loop integral can be read off directly from the 
underlying Feynman graph. The first part of this course closes with 6, which shows how finite 
results are obtained within perturbation theory. The remaining section are more mathematical 
in nature: Section 7 states a general theorem which relates Feynman integrals to periods. Shuf- 
fle algebras and multiple polylogarithms are treated in section 8 and section 9, respectively. In 
sect. 10 we discuss how multiple polylogarithms emerge in the calculation of Feynman integrals. 
Finally, section 1 1 provides a summary. 

2 Basics of perturbative quantum field theory 

Elementary particle physics is described by quantum field theory. To begin with let us start with 
a single field ^{x) . Important concepts in quantum field theory are the Lagrangian, the action and 
the generating functional. If (|)(x) is a scalar field, a typical Lagrangian is 

L = l(a,(^(x))(a^(^w)-lm>w2 + lx^W4. (2) 

The quantity m is interpreted as the mass of the particle described by the field ^{x), the quantity 
X describes the strength of the interactions among the particles. Integrating the Lagrangian over 
Minkowski space yields the action: 

5[^] = fd^xL{(^). (3) 

The action is a functional of the field (|). In order to arrive at the generating functional we in- 
troduce an auxiliary field J{x), called the source field, and integrate over all field configurations 



Z[J] = 5V; /"2)(^e'(^['^]+/'^'^-/W'^W). (4) 

The integral over all field configurations is an infinite-dimensional integral. It is called a path 
integral. The prefactor 9{, is chosen such that Z[0] = 1. The n-point Green function is given by 

With the help of functional derivatives this can be expressed as 

5«Z[7] 



{0\T{<^{xi)..Mxn))\0) = {-if 



hJ{xi)...hJ{Xn) 



(6) 



We are in particular interested in connected Green functions. These are obtained from a func- 
tional W[J], which is related to Z[J] by 

Z[J] = e^^^\ (7) 

The connected Green functions are then given by 



G„(xi,...,x„) = (-/) 



..^„-l 5"iy[7] 



hJ{xi)...hJ{Xn) 



(8) 
y=o 



It is convenient to go from position space to momentum space by a Fourier transformation. We 
define the Green functions in momentum space by 

Gn{xu...,xn) = l^...^e-^^P^^^{2n)'d{p, + ...+p,)Gn{pu....Pn). (9) 

Note that the Fourier transform G„ is defined by explicitly factoring out the 5-function 5(pi + 
... + Pn) and a factor {2n)^. We denote the two-point function in momentum space by G2{p). In 
this case we have to specify only one momentum, since the momentum flowing into the Green 
function on one side has to be equal to the momentum flowing out of the Green function on the 
other side due to the presence of the 5-function in eq. (9) . We now are in a position to define the 
scattering amplitude: In momentum space the scattering amplitude with n external particles is 
given by the connected n-point Green function multiplied by the inverse two-point function for 
each external particle: 

Jln{PU---,Pn) = GlipiY ...GlipnY Gn{pu ..., Pn)- (10) 

The scattering amplitude enters directly the calculation of a physical observable. Let us first 
consider the scattering process of two incoming particles with four-momenta p[ and pj and 
{n — 2) outgoing particles with four- momenta pi to Pn-i- Let us assume that we are interested 



in an observable 0{pi, ....Pn-i) which depends on the momenta of the outgoing particles. In 
general the observable depends on the experimental set-up and can be an arbitrary complicated 
function of the four-momenta. In the simplest case this function is just a constant equal to one, 
corresponding to the situation where we count every event with (n — 2) particles in the final 
state. In more realistic situations one takes for example into account that it is not possible to 
detect particles close to the beam pipe. The function O would then be zero in these regions 
of phase space. Furthermore any experiment has a finite resolution. Therefore it will not be 
possible to detect particles which are very soft or which are very close in angle to other particles. 
We will therefore sum over the number of final state particles. In order to obtain finite results 
within perturbation theory we have to require that in the case where one or more particles become 
unresolved the value of the observable has a continuous limit agreeing with the value of the 
observable for a configuration where the unresolved particles have been merged into "hard" (or 
resolved) particles. Observables having this property are called infrared-safe observables. The 
expectation value for the observable O is given by 

(<^) = 0/ / I M2 E/^^n-2Q(Pl,-,Pn-2)|-^«P, (H) 

APl+P2)nJ 

where 1/2/ {p[ +^2)^ i^ ^ normalisation factor taking into account the incoming flux. The phase 
space measure is given by 



n\ fj^ {InyiEi 



i=i 



Ei is the energy of particle /: 



?2 I ,^2 



Ei ^ \^pf + mf (13) 



We see that the expectation value of O is given by the phase space integral over the observable, 
weighted by the norm squared of the scattering amplitude. As the integrand can be a rather 
complicated function, the phase space integral is usually performed numerically by Monte Carlo 
integration. 

Let us now look towards a more realistic theory. As an example I will take quantum chromo- 
dynamics (QCD), which describes the strong force and which is formulated in terms of quarks 
and gluons. Quarks and gluons are collectively called partons. There are a few modifications to 
eq. (11). The master formula reads now 

(O) = £ fdxifaiXi) [dX2fb{x2) (14) 

a,b-' -^ 

2^n,(l)n,(2)n,(l)n,(2) ^J spmsT^lour 

The partons have internal degrees of freedom, given by the spin and the colour of the partons. In 
squaring the amplitude we sum over these degrees of freedom. For the particles in the initial state 



we would like to average over these degrees of freedom. This is done by dividing by the factors 
ns{i) and nc{i), giving the number of spin degrees of freedom (2 for quarks and gluons) and the 
number of colour degrees of freedom (3 for quarks, 8 for gluons). The second modification is 
due to the fact that the particles brought into collision are not partons, but composite particles 
like protons. At high energies the constituents of the protons interact and we have to include a 
function fa{x) giving us the probability of finding a parton a with momentum fraction x of the 
original proton momentum inside the proton, s is the centre-of-mass energy squared of the two 
partons entering the hard interaction. In addition there is a small change in eq. (12). The quantity 
(n!) is replaced by (H^j 0' where rij is the number of times a parton of type j occurs in the final 
state. 

As before, the scattering amplitude Rn can be calculated once the Lagrangian of the theory 
has been specified. For QCD the Lagrange density reads: 



^ quarks q 

with 

F-,{x) = d,A'^,{x) -avA«(..) +^r*^Aj(..)A^, D, = d,- igT-A-{x). (16) 

The gluon field is denoted by A"{x), the quark fields are denoted by \\fq{x). The sum is over all 
quark flavours. The masses of the quarks are denoted by niq. There is a summation over the 
colour indices of the quarks, which is not shown explicitly. The variable g gives the strength of 
the strong coupling. The generators of the group SU{3) are denoted by T" and satisfy 



rpa rpb 



if^'T'. (17) 



The quantity F^ is called the field strength, the quantity D^ is called the covariant derivative. 
The variable ^ is called the gauge-fixing parameter. Gauge-invariant quantities like scattering 
amplitudes are independent of this parameter. L^v stands for the Faddeev-Popov term, which 
arises through the gauge-fixing procedure and which is only relevant for loop amplitudes. 

Unfortunately it is not possible to calculate from this Lagrangian the scattering amplitude Jl„ 
exactly. The best what can be done is to expand the scattering amplitude in the small parameter g 
and to calculate the first few terms. The amplitude ^n with n external partons has the perturbative 
expansion 

^n = /-^f^r+/^i^)+AP+/^F + ...V (18) 



In principle we could now calculate every term in this expansion by taking the functional deriva- 
tives according to eq. (8). This is rather tedious and there is a short-cut to arrive at the same 
result, which is based on Feynman graphs. The recipe for the computation of Jl„ is as follows: 
Draw first all Feynman diagrams with the given number of external particles and / loops. Then 
translate each graph into a mathematical formula with the help of the Feynman rules, ^n is then 
given as the sum of all these terms. 



In order to derive the Feynman rules from the Lagrangian one proceeds as follows: One first 
separates the Lagrangian into a part which is bilinear in the fields, and a part where each term 
contains three or more fields. (A "normal" Lagrangian does not have parts with just one or zero 
fields.) From the part bilinear in the fields one derives the propagators, while the terms with 
three or more fields give rise to vertices. As an example we consider the gluonic part of the QCD 
Lagrange density: 

^QCD = 



Ik^x) 




A^x 


-gr'^ (a^A^(x)) A^^(^)A^-^(x) - ^-g^r'r^'. 


T"-^ quarks 


+ iFP. 





(19) 

Within perturbation theory we always assume that all fields fall off rapidly enough at infinity. 
Therefore we can ignore boundary terms within partial integrations. The expression in the first 
line is bilinear in the fields. The terms in the square bracket in this line define an operator 

pMVab^^^ = dpdPg^''d"''-(l-^^d^d''d"'\ (20) 

For the propagator we are interested in the inverse of this operator 

p,aac^^^^p-iy^^^_y^ = g%d"'d\x-y). (21) 



Working in momentum space we are more specifically interested in the Fourier transform of the 
inverse of this operator: 

^iv( 
the gluon propagator one finds the Feynman rule 



The Feynman rule for the propagator is then given by (P ^)uv(^) times the imaginary unit. For 



^'^vimm.^'^ = ^(-^,v + (l-^)^)6"^ (23) 

To derive the Feynman rules for the vertices we look as an example at the first term in the second 
line of eq. (19): 

^ggs = -gr''{d,AUx))A'^{x)A^'ix). (24) 

This term contains three gluon fields and will give rise to the three-gluon vertex. We rewrite this 
term as follows: 

Lggg = ld%d%d''x3a'''''^"'^{x,xuX2,X3)A;{xi)A'^{x2)A{{x3), (25) 



where 



a"^'^''\x,x,,X2,x^) = gr''-g^^^{d\h\x-x,))h\x-X2)h\x-x,). (26) 



Again we are interested in the Fourier transform of this expression: 

d'^ki d'^k2 d^k^i _■ 



a 



abc fjvX 



(a:,Xi,X2,JC3) 



/ 



(271)4(271)4(271)4 



iki (xi —x) — ik2 {x2 —x)— ik^ {xj —x) f^abc /jV A. 



Working this out we find 



ce^'^^^\hMM] 



-gf^'g^'^iki 



"{ki.kiM)- 



(27) 



The Feynman rule for the vertex is then given by the sum over all permutations of identical 
particles of the function a multiplied by the imaginary unit i. (In the case of identical fermions 
there would be in addition a minus sign for every odd permutation of the fermions.) We thus 
obtain the Feynman rule for the three-gluon vertex: 




gf 



■abc 



r (k\ - k\) + /^ (^ - ^) + /^ {k\ - kl) 



(28) 



Note that there is momentum conservation at each vertex, for the three-gluon vertex this implies 

ki+k2+k^ = 0. (29) 

Following the procedures outlined above we can derive the Feynman rules for all propagators 
and vertices of the theory. If an external particle carries spin, we have to associate a factor, which 
describes the polarisation of the corresponding particle when we translate a Feynman diagram 
into a formula. Thus, there is a polarisation vector £^(/?) for each external gauge boson and a 
spinor u{p), u{p), v{p) or v{p) for each external fermion. 

Furthermore there are a few additional rules: First of all, there is an integration 



d^k 

(2^ 



(30) 



for each internal momentum which is not constrained by momentum conservation. Such an 
integration is called a "loop integration" and the number of independent loop integrations in a 
diagram is called the loop number of the diagram. Secondly, each closed fermion loop gets an 
extra factor of (—1). Finally, each diagram gets multiplied by a symmetry factor 1/S, where 
S is the order of the permutation group of the internal lines and vertices leaving the diagram 
unchanged when the external lines are fixed. 



Let us finish this section by listing the remaining Feynman rules for QCD. The quark and the 
ghost propagators are given by 



, 1^ + m 



5, 



F-m2"^'' 



The Feynman rules for the four-gluon vertex, the quark-gluon vertex and the ghost-gluon vertex 
are 




-ig" 


■fabefecd^^,X^vp_g,Pg 


^j-adefetc(^g,.gXp_g,XgVpJ 


igfTfi. 



^^^ , = -gf'^'k^- (32) 

*' 

q,c' 

The Feynman rules for the electro-weak sector of the Standard Model are similar, but too numer- 
ous to list them explicitly here. 

Having stated the Feynman rules, let us look at some examples. We have seen that for a 
given process with a specified set of external particles the scattering amplitude is given as the 
sum of all Feynman diagrams with this set of external particles. We can order the diagrams 
by the powers of the coupling factors. In QCD we obtain for each three-particle vertex one 
power of g, while the four-gluon vertex contributes two powers of g. The leading order result 
for the scattering amplitude is obtained by taking only the diagrams with the minimal number 
of coupling factors g into account. These are diagrams which have no closed loops. There are 
no conceptual difficulties in evaluating these diagrams. However going beyond the leading order 
in perturbation theory, loop diagrams appear which involve integrations over the loop momenta. 
These diagrams are more difficult to evaluate and I will discuss them in more detail. Fig. 1 shows 
a Feynman diagram contributing to the one-loop corrections for the process e+e — )■ qgq. At high 
energies we can ignore the masses of the electron and the light quarks. From the Feynman rules 
one obtains for this diagram: 

-eVCFr,';v>4)A(;.5)^/^^"bi)?'b2)^Yv^Y,§Y^(P3). (33) 

Pj23 J l^'^J K'2 Pl2 '^l % 




P4 / \ P3 

Figure 1: A one-loop Feynman diagram contributing to the process e^e~ — )• qgq. 

Here, pn = Pi +P2, Pm = Pi +P2 + P3, h = h -pn, ^3 = ^2 -P3- Further ^(^2) = YTe^(P2), 
where e^(p2) is the polarisation vector of the outgoing gluon. All external momenta are assumed 
to be massless: pj = for i = 1..5. We can reorganise this formula into a part, which depends 
on the loop integration and a part, which does not. The loop integral to be calculated reads: 

while the remainder, which is independent of the loop integration is given by 

-e'g^CFTfiv{p^)fu{p5)^^u{p,)4{p2ynlvl^l^r{<^fv{p^). (35) 

PiiifPn 

The loop integral in eq. (34) contains in the denominator three propagator factors and in the 
numerator two factors of the loop momentum. We call a loop integral, in which the loop momen- 
tum occurs also in the numerator a "tensor integral". A loop integral, in which the numerator is 
independent of the loop momentum is called a "scalar integral". The scalar integral associated to 
eq. (34) reads 

d'ki 1 

{2TiYk\kikY ^ ^ 

It is always possible to reduce tensor integrals to scalar integrals [48,49]. The calculation of 
integrals like the one in eq. (36) is the main topic of these lectures. More information on the 
basics of perturbation theory and quantum field theory can be found in one of the many textbooks 
on quantum field theory, like for example in refs. [50, 51]. 

3 Dimensional regularisation 

Before we start with the actual calculation of loop integrals, I should mention one complication: 
Loop integrals are often divergent ! Let us first look at the simple example of a scalar two-point 
one-loop integral with zero external momentum: 



JU 1 1 /• ..-, 1 1 rdx 




1 /"^t2 1 _ 1 fdx 



(271)4 (/t2)2 (471)2; k^ (471) 
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This integral diverges at fc^ — )> oo as well as at k^ — )> 0. The former divergence is called ultraviolet 
divergence, the later is called infrared divergence. Any quantity, which is given by a divergent 
integral, is of course an ill-defined quantity. Therefore the first step is to make these integrals 
well-defined by introducing a regulator. There are several possibilities how this can be done, but 
the method of dimensional regularisation [52-54] has almost become a standard, as the calcula- 
tions in this regularisation scheme turn out to be the simplest. Within dimensional regularisation 
one replaces the four-dimensional integral over the loop momentum by an Z)-dimensional inte- 
gral, where D is now an additional parameter, which can be a non-integer or even a complex 
number. We consider the result of the integration as a function of D and we are interested in 
the behaviour of this function as D approaches 4. The D-dimensional integration still fulfils the 
standard laws for integration, like linearity, translation invariance and scaling behaviour [55,56]. 
If / and g are two functions, and if a and b are two constants, linearity states that 

fd^k{af{k)+bg{k)) = a fd^kf{k)+b fd^kg{k). (38) 

Translation invariance requires that 

fd^kf{k + p) = fd^kfik). (39) 

for any vector p. 

The scaling law states that 

fd'^kfiXk) = X-^ fd^kfik). (40) 

The D-dimensional integral has also a rotation invariance: 

fd^kf{Ak) = fd^kfik), (41) 

where A is an element of the Lorentz group SO{l,D— 1) of the D-dimensional vector-space. 
Here we assumed that the D-dimensional vector-space has the metric diag(+l, — 1, — 1, — 1, ...). 
The integral measure is normalised such that it agrees with the result for the integration of a 
Gaussian function for all integer values D: 

/ D 

d'^kcxp{ak^) = ii-Y ■ (42) 

We will further assume that we can always decompose any vector into a 4-dimensional part and 
a (D — 4) -dimensional part 

and that the 4-dimensional and (Z) — 4) -dimensional subspaces are orthogonal to each other: 

k(^4)-k(D-4) = 0. (44) 

11 



If D is an integer greater than 4, this is obvious. We postulate that these relations are true for any 
value of D. One can think of the underlying vector-space as a space of infinite dimension, where 
the integral measure mimics the one in D dimensions. 

In practise we will always arrange things such that every function we integrate over D dimen- 
sions is rotational invariant, e.g. is a function of k^. In this case the integration over the (D — 1) 
angles is trivial and can be expressed in a closed form as a function of D. Let us assume that we 
have an integral, which has a UV-divergence, but no IR-divergences. Let us further assume that 
this integral would diverge logarithmically, if we would use a cut-off regularisation instead of 
dimensional regularisation. It turns out that this integral will be convergent if the real part of D 
is smaller than 4. Therefore we may compute this integral under the assumption that Re(D) < 4 
and we will obtain as a result a function of D. This function can be analytically continued to the 
whole complex plane. We are mainly interested in what happens close to the point D = 4. For 
an ultraviolet divergent one-loop integral we will find that the analytically continued result will 
exhibit a pole at D = 4. It should be mentioned that there are also integrals which are quadrati- 
cally divergent, if a cut-off regulator is used. In this case we can repeat the argumentation above 
with the replacement Re(D) < 2. 

Similarly, we can consider an IR-divergent integral, which has no UV-divergence. This in- 
tegral will be convergent if Re(D) > 4. Again, we can compute the integral in this domain and 
continue the result to D = 4. Here we find that each IR-divergent loop integral can lead to a 
double pole at D = 4. 

We will use dimensional regularisation to regulate both the ultraviolet and infrared diver- 
gences. The attentative reader may ask how this goes together, as we argued above that UV- 
divergences require Re(D) < 4 or even Re(Z)) < 2, whereas IR-divergences are regulated by 
Re(D) > 4. Suppose for the moment that we use dimensional regularisation just for the UV- 
divergences and that we use a second regulator for the IR-divergences. For the IR-divergences 
we could keep all external momenta off-shell, or introduce small masses for all massless parti- 
cles or even raise the original propagators to some power v. The exact implementation of this 
regulator is not important, as long as the IR-divergences are screened by this procedure. We 
then perform the loop integration in the domain where the integral is UV-convergent. We ob- 
tain a result, which we can analytically continue to the whole complex D-plane, in particular to 
Re(Z)) > 4. There we can remove the additional regulator and the IR-divergences are now regu- 
lated by dimensional regularisation. Then the infrared divergences will also show up as poles at 
D = 4. 

There is one more item which needs to be discussed in the context of dimensional regulari- 
sation: Let us look again at the example in eqs. (33) to (35). We separated the loop integral from 
the remainder in eq. (35), which is independent of the loop integration. In this remainder the 
following string of Dirac matrices occurs: 

YvYpY^YaY". (45) 

If we anti-commute the first Dirac matrix, we can achieve that the two Dirac matrices with index 
V are next to each other: 

YvY"- (46) 

12 



In four dimensions this equals 4 times the unit matrix. What is the value within dimensional 
regularisation ? The answer depends on how we treat the Dirac algebra. Does the Dirac algebra 
remain in four dimensions or do we also continue the Dirac algebra to D dimensions ? There 
are several schemes on the market which treat this issue differently. To discuss these schemes it 
is best to look how they treat the momenta and the polarisation vectors of observed and unob- 
served particles. Unobserved particles are particles circulating inside loops or emitted particles 
not resolved within a given detector resolution. The most commonly used schemes are the con- 
ventional dimensional regularisation scheme (CDR) [56], where all momenta and all polarisation 
vectors are taken to be in D dimensions, the 't Hooft-Veltman scheme (HV) [52, 57], where the 
momenta and the helicities of the unobserved particles are D dimensional, whereas the momenta 
and the helicities of the observed particles are 4 dimensional, and the four-dimensional helicity 
scheme (FD) [58-60], where all polarisation vectors are kept in four dimensions, as well as the 
momenta of the observed particles. Only the momenta of the unobserved particles are continued 
to D dimensions. 

The conventional scheme is mostly used for an analytical calculation of the interference of 
a one-loop amplitude with the Bom amplitude by using polarisation sums corresponding to D 
dimensions. For the calculation of one-loop helicity amplitudes the 't Hooft-Veltman scheme and 
the four-dimensional helicity scheme are possible choices. All schemes have in common, that the 
propagators appearing in the denominator of the loop-integrals are continued to D dimensions. 
They differ how they treat the algebraic part in the numerator. In the 't Hooft-Veltman scheme the 
algebraic part is treated in D dimensions, whereas in the FD scheme the algebraic part is treated 
in four dimensions. It is possible to relate results obtained in one scheme to another scheme, 
using simple and universal transition formulae [61-63]. Therefore, if we return to the example 
above, we have 

V _ r D-l, in the CDR and HV scheme, 
'^vT - I 4-1, in the FD scheme. ^ ' 

To summarise we are interested into loop integrals regulated by dimensional regularisation. As 
a result we seek the Laurent expansion around D = A. It is common practise to parametrise the 
deviation of D from 4 by 

D = 4-2£. (48) 

Divergent loop integrals will therefore have poles in 1/e. In an Z-loop integral ultraviolet diver- 
gences will lead to poles 1 /e' at the worst, whereas infrared divergences can lead to poles up to 

l/£2'. 

4 Loop integration in D dimensions 

In this section I will discuss how to perform the D-dimensional loop integrals. It would be more 
correct to say that we exchange them for some parameter integrals. As an example we take the 
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one-loop integral of eq. (36): 



I = f^:^ I (49) 



The integration is now in D dimensions. In eq. (49) there are some overall factors, which I 
inserted for convenience: The integral measure is now d^k/{i%^'^) instead of d^k/{2%)^, and 
each propagator is multiplied by (—1). The reason for doing this is that the final result will be 
simpler. 

As already discussed above, the only functions we really want to integrate over D dimensions 
are the ones which depend on the loop momentum only through k^. The integrand in eq. (49) 
is not yet in such a form. To bring the integrand into this form, we first convert the product 
of propagators into a sum. We can do this with the Feynman parameter technique. In its full 
generality it is also applicable to cases, where each factor in the denominator is raised to some 
power V. The formula reads: 



tlAy. = ^^^j in^^^-^ri / /^\^v , v = £v,. (50) 

The proof of this formula can be found in many text books and is not repeated here. r{x) is 
Euler's Gamma function, 5(x) denotes Dirac's delta function. The price we have to pay for 
converting the product into a sum are (n — l) additional integrations. Let us look at the example 
from eq. (36): 

1 1 1 

2 /"j- /"j- /"j- 5(l-JCi-JC2-JC3) 




/rf.,/</../</x3 -^'-^^_-2 (51) 



In the next step we complete the square and shift the loop momentum, such that the integrand 
becomes a function of k^. With k2 = k\— pn and % = A:2 — ps we have 

~Xxk\-X2kl-Xj,kl = -{ki-X2Pl2-X3Pi23) -■«1-«2*12--«1-«3>^123, (52) 

where 512 = {pi + pi)^ and 5123 = {pi + P2 + Ps)^- We can now define 

k[ = ki-X2Pn-X3Pi23 (53) 

and using translational invariance our loop integral becomes 

(-^1 -X1X2S12-X1X3S123J 

/ 2 7 

The integrand is now a function of k[ , which we can relabel as k . 
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Im ^0 




Re ^0 



Figure 2: Integration contour for the Wick rotation. The little circles along the real axis exclude 
the poles. 

Having succeeded to rewrite the integrand as a function of k^, we then perform a Wick rota- 
tion, which transforms Minkowski space into an Euclidean space. Remember, that k^ written out 
in components in D-dimensional Minkowski space reads 



kl-kj-kj-kj-... 



(55) 



(Here kj denotes the j-th component of the vector k, in contrast to the previous notation, where 
we used the subscript to label different vectors kj. It should be clear from the context what 
is meant.) Furthermore, when integrating over kg, we encounter poles which are avoided by 
Feynman's ?'5-prescription. In the complex ^o-plane we consider the integration contour shown 
in fig. 2. Since the contour does not enclose any poles, the integral along the complete contour 
is zero: 



jxikofiko) = 0. 



(56) 



If the quarter-circles at infinity give a vanishing contribution (it can be shown that this is the case) 
we obtain 



dkof{ko) 



dkofiko) 



(57) 



We now make the following change of variables: 



Kj, forl<j<D-l. 



(58) 



As a consequence we have 



k^ = -K^, d'^k = id^K, 



(59) 
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where K is now given with Euclidean signature: 

K^ = kI + kI+kI+kI + ... (60) 

Combining eq. (57) with eq. (58) we obtain for the integration of a function f{k^) over D dimen- 
sions 

whenever there are no poles inside the contour of fig. 2 and the arcs at infinity give a vanishing 
contribution. The integral on the r.h.s. is now over D-dimensional Euclidean space. Eq. (61) 
justifies our conventions, to introduce a factor i in the denominator and a minus sign for each 
propagator in eq. (49). These conventions are just such that after Wick rotation we have simple 
formulae. 

We now have an integral over D-dimensional Euclidean space, where the integrand depends 
only on K^. It is therefore natural to introduce spherical coordinates. In D dimensions they are 
given by 

^0 = ^cosGi, 

^1 = ^sin0icos02, 

Kd-2 = ^sin0i...sin0D 2COS0D 1, 

Kd I = ^sin0i...sin0D_2sin0£)_i. (62) 

In D dimensions we have one radial variable K, D — 2 polar angles ^ (with 1 < j < D — 2) and 
one azimuthal angle 0o-i. The measure becomes 



d'^K = K^^dKdQD, J^o= n^in^^^^'Q'^Q'- (63) 



271^/2 



D-1 

n 

1=1 

Integration over the angles yields 

IT 71 231 

fd^D = f dBisin'^-^Bi... f dQD-2sinQD-2 f dBD-i = ^r^. (64) 

^^' 

Note that the integration on the l.h.s of eq. (64) is defined for any natural number D, whereas the 
result on the r.h.s is an analytic function of D, which can be continued to any complex value. 

It is now the appropriate place to say a few words on Euler's Gamma function. The Gamma 
function is defined for Re(jc) > by 

r{x) = / e-'f-'^dt. (65) 
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It fulfils the functional equation 

r(;c+l) = xr{x). (66) 

For positive integers n it takes the values 

r(n + l) = n! = l-2-3-...-n. (67) 

For integers n we have the reflection identity 

^x-n) r(i-.y) 

r{x) r{\ —x+n) 

The Gamma function r{x) has poles located on the negative real axis atx = 0,-1,-2, .... Quite 
often we will need the expansion around these poles. This can be obtained from the expansion 
around x= I and the functional equation. The expansion around £ = 1 reads 

r(l+8) = expf-Y££+£^^^„£"J, (69) 

where Je is Euler's constant 

JE = lim I y^-lnn) =0.5772156649... (70) 



and C,n is given by 



OO 1 

^n = L-n- (71) 

For example we obtain for the Laurent expansion around £ = 

r(e) = ^-Y£ + 0(£). (72) 

We are now in a position to perform the integration over the loop momentum. Let us discuss 
again the example from eq. (54). After Wick rotation we have 

d^ki 1 _9 M^^ /"jS ^l-Xi-X2-X3) 



mDlH-k\){-kl){-kl) J n^J {K^-xiX2Sn-xix,sn3 



Introducing spherical coordinates and performing the angular integration this becomes 



(73) 



h{\-xi-X2~X3){K^y 



D-l 



2 /■ 9 /" o b[\-x\-X2~x-i)\ti'-] "■ 
I = -j^ / dK^ / d^x- '^ ^ 3 . (74) 



^(t)^ ^ {K^-xiX2Sn-xiX2,sn2,y 
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For the radial integration we have after the substitution t = K^/{—x\X2S\2 — ^i^3i'i23) 



D 2 oo £,_2 



dK T = {-X]_X2Si2-XiX'iSi2'i)^ dt 



\3' 



(75) 



{K^-XiX2Sl2-XiX3Si23) { (1+0 



The remaining integral is a standard integral and yields 

dt ^ = \^L —■ (76) 

(1+0' r(3) 

Putting everything together and setting D = 4 — 2£ we obtain 

/ = r(l+.)/A8(l-x,-..-X3)xr-(-X2.,2-X3.„3)-'-. (77) 

Therefore we succeeded in performing the integration over the loop momentum k at the expense 
of introducing a two-fold integral over the Feynman parameters. We will learn techniques how 
to perform the Feynman parameter integrals later in these lectures. Let me however already state 
the final result: 



1 



^123-^12 



( - - Y£ - In (-5123) J ln.)c - - In^^jc 



+ 0(£), x = ^^. (78) 

-5123 



The result has been expanded in the regularisation parameter £ up to the order 0(£). We see 
that the result has a term proportional to l/£. Poles in £ in the final (regularised) result reflect 
the original divergences in the unregularised integral. In this example the pole corresponds to a 
coUinear singularity. 

5 Multi-loop integrals 

As the steps discussed in the previous section always occur in any loop integration we can com- 
bine them into a master formula. Let us consider a scalar Feynman graph G with m external lines 
and n internal lines. We denote by Iq the associated scalar /-loop integral. For each internal line 
j the corresponding propagator in the integrand can be raised to an integer power Vj. Therefore 
the integral will depend also on the numbers Vi,...,V„. 

J r=\ 1112 j^i {-qj + nij) J 

The independent loop momenta are labelled k\, ...,ki. The momenta flowing through the propa- 
gators are then given as a linear combination of the external momenta p and the loop momenta k 
with coefficients — 1 , or 1 : 

I m 

qt = Y.^ijkj + lL^ijPh hj^Oije {-1,0,1} . (80) 

7=1 7=1 
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Pi 



PA 



Figure 3: The "double box"-graph: A two-loop Feynman diagram with four external lines and 
seven internal lines. The momenta flowing out along the external lines are labelled p\, ..., P4, the 
momenta flowing through the internal lines are labelled qi, ..., qq. 

We can repeat for each loop integration the steps of the previous section. Doing so, we arrive at 
the following Feynman parameter integral: 



r(v-ZD/2) /• ,„ ^,, V ^/^T^ v,-i\W 



V-(/+l)D/2 



nr(v,) ^/>o 

7=1 



i=\ \j=\ 



jV-lD/l 



(81) 



The functions 11 and iT depend on the Feynman parameters Xj. If one expresses 



;=i 



£ £ KMrsh + £ 2A:, ■ 2, + 7, 

r=\s=\ r=\ 



(82) 



where M is a Z x / matrix with scalar entries and 2 is a /-vector with four-vectors as entries, one 
obtains 



U = det(M), J = det(M) {J + QM-^Q) . 



(83) 



As an example let us look at the two-loop double box graph of fig. (3). In fig. 3 there are two 
independent loop momenta. We can choose them to be ki = q^ and ^2 = <?6- Then all other 
internal momenta are expressed in terms ofki, ki and the external momenta />i, ..., ;?4: 



qi=ki-pi, q2 = ki-pi-p2, <?4 = ^l+^2, 

q5 = k2-p3-P4, qi^ki-pA- 



(84) 



We will consider the case 



p\ = Q, pI = 0, pI=0, pI = 0, 

mi = m2 = niT, = m^ = m^ = mg = m-j = 0. 



(85) 



We define 



s = {pi+Pif = iP3+P4)^, t = {pj+p-if- = ipi+P4) 



(86) 
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We have 



E-^j(-^?) 



(^1 +X2 + X3 +X4) k\ - Ix^ki ■ k2 - {x4 +X5 +xe + jc?) ^2 



(87) 



+2 [Xipi +X2 {pi +P2)] -ki+l [X5 {p3 + Pa) +XtP4] ■ ^2 - (-^2 +-^5) *• 
In comparing with eq. (82) we find 



M 

Q 
J 



Xi+X2+X2,+X4 X4 

X4 X4 + X5 + X(, + X-] 

XiPi+X2{pi+P2) 
X5ip3+P4)+X7P4 

{X2+X5){-S). 



(88) 



Plugging this into eq. (83) we obtain the graph polynomials as 

Zl = {xl+X2+X3){x5+X(,+X^)+X4{xi+X2+X3+X5+Xe+X■^), 

7 = [jC2JC3(-^4+-^5+-^6+-^7)+-«5-^6(-^l +-^2+-^3+-«4) +-^2-«4-^6+-^3-«4-«5] (-■^) 
+X[X4X'] {—t) . 



(89) 



There are several other ways how the two polynomials U and 7 can be obtained [64]. Let me 
mention one method, where the two polynomials can be read off directly from the topology of 
the graph G. We consider first connected tree graphs T, which are obtained from the graph 
G by cutting / lines. The set of all such trees (or 1 -trees) is denoted by Ti. The Feynman 
parameters corresponding to the cut lines define a monomial of degree I. 11 is the sum over all 
such monomials. Cutting one more line of a 1-tree leads to two disconnected trees {Ti,T2), or 
a 2-tree. T2 is the set of all such pairs. The cut lines define monomials of degree / + 1. Each 
2-tree of a graph corresponds to a cut defined by cutting the lines which connected the two now 
disconnected trees in the original graph. The square of the sum of momenta through the cut lines 
of one of the two disconnected trees Ti or T2 defines a Lorentz invariant 



HtuTi) 



kM{TuT2) 



(90) 



The function ^q is the sum over all such monomials times minus the corresponding invariant. 
The function J is then given by ^0 plus an additional piece involving the internal masses nij. In 
summary, the functions U and J are obtained from the graph as follows: 



TeTi jiT 

i- , n 

{Ti,T2)eT2 mTi,T2) J=l 



(91) 
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6 How to obtain finite results 

We have already seen in eq. (78) that the final result of a regularised Feynman integral may 
contain poles in the regularisation parameter £. These poles reflect the original ultraviolet and 
infrared singularities of the unregularised integral. What shall we do with these poles ? The 
answer has to come from physics and we distinguish again the case of UV-divergences and IR- 
divergences. The UV-divergences are removed through renormalisation. Ultraviolet divergences 
are absorbed into a redefinition of the parameters. As an example we consider the renormalisation 
of the coupling: 

g = Zg g, . (92) 

divergent divergent finite 

The renormalisation constant Zg absorbs the divergent part. However Zg is not unique: One may 
always shift a finite piece from gy to Zg or vice versa. Different choices for Zg correspond to dif- 
ferent renormalisation schemes. Two different renormalisation schemes are always connected by 
a finite renormalisation. Note that different renormalisation schemes give numerically different 
answers. Therefore one always has to specify the renormalisation scheme. Some popular renor- 
malisation schemes are the on-shell scheme, where the renormalisation constants are defined by 
conditions at a scale where the particles are on-shell. A second widely used scheme is modified 
minimal subtraction. In this scheme one always absorbs the combination 

A = --Y£+ln47r (93) 

into the renormalisation constants. One proceeds similar with all other quantities appearing in 
the original Lagrangian. For example: 

A"^ = v^A;^ ,, \\f^ = VZ^\\fq,r, g = Zggr, m = Zmntr, ^ = Z^^,. (94) 

The fact that square roots appear for the field renormalisation is just convention. Let us look a 
little bit closer into the coupling renormalisation within dimensional regularisation and the MS- 
renormalisation scheme. Within dimensional regularisation the renormalised coupling gr is a 
dimensionfuU quantity. We define a dimensionless quantity g^ by 

gr = gRlU^ (95) 

where fj is an arbitrary mass scale. From a one-loop calculation one obtains 

Zs = l-\^ojl§^A+0{gt), Po = yA^c - ^A^/. (96) 

Nc is the number of colours and Nf the number of light quarks. The quantity gR will depend on 
the arbitrary scale /u. To derive this dependence one first notes that the unrenormalised coupling 
constant g is of course independent of /j: 

^g = (97) 
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Substituting g = Zgf/gR into this equation one obtains 



From eq. (96) one obtains 



V' = -^^'^-{^s'^T/']^^- 



'^'4/^ - p«(4+^(^^) 



(98) 



(99) 



Instead of ^^ one often uses the quantity a^ = g^/{4n). Going to D = 4 one arrives at 

i^ = .^(^)\„(^V_ (100) 



'""'•- Msj +^4./ 



df?4it 



This differential equation gives the dependence of a^ on the renormalisation scale fu. At leading 
order the solution is given by 



«.v(A') 
471 



1 



Poln^ 



-,2\ ' 



(101) 



where A is an integration constant. The quantity A is called the QCD scale parameter. For QCD 
Po is positive and CLs{i-i) decreases with larger /j. This property is called asymptotic freedom: 
The coupling becomes smaller at high energies. In QED |3o has the opposite sign and the fine- 
structure constant a(/j) increases with larger /j. The electromagnetic coupling becomes weaker 
when we go to smaller energies. 

Let us now look at the infrared divergences: We first note that any detector has a finite 
resolution. Therefore two particles which are sufficiently close to each other in phase space will 
be detected as one particle. Now let us look again at eqs. (1 1) and (18). The next-to-leading order 
term will receive contributions from the interference term of the one-loop amplitude Jin with the 
leading-order amplitude ^„ , both with {n — 2) final state particles. This contribution is of order 
^2m-2 q^ j-j^g same order is the square of the leading-order amplitude ■^lj_i with (n — l) final 
state particles. This contribution we have to take into account whenever our detector resolves 
only n particles. It turns out that the phase space integration over the regions where one or more 
particles become unresolved is also divergent, and, when performed in D dimensions, leads to 
poles with the opposite sign as the one encountered in the loop amplitudes. Therefore the sum of 
the two contributions is finite. The Kinoshita-Lee-Nauenberg theorem [65,66] guarantees that all 
infrared divergences cancel, when summed over all degenerate physical states. As an example 
we consider the NLO corrections to y* — )> 2 jets. The interference term of the one-loop amplitude 
with the Born amplitude is given by 



2Re Jl 



(0)* 



-^o 



OC.S 

71 



Cf 



1 



3 . 7 2 
:^-4 + --7r^ 
2£ 12 



R 



(0) 



+ 0(£). 



(102) 
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Se = (47r)^e^^^^ is the typical phase-space volume factor in D = 4 — 2£ dimensions. For sim- 
plicity we have set the renormalisation scale /u equal to the centre-of-mass energy squared s. The 
square of the Born amplitude is given by 



Jl 



(0) 



= 167rA/ca(l -£)i'. 



(103) 



This is independent of the final state momenta and the integration over the phase space can be 
written as 



/ d^2 h 



d<^2 (2 Re -^f^*-^3^' 



a. 



'-Cr 



1 



3_ 
2e 



, 7 r 



d(^2 



SI 



(0) 



+ 0(£). (104) 



71 V £ 
The real corrections are given by the leading order matrix element for y* — )■ qgq and read 



SI 



(0) 



n^n^aasCFNcil-t) 



^---- + (1 

X\X2 X\ X2 



.e)^ + (l-e)^-2e 

X\ X2 



(105) 



where x\ = 512/5123, X2 = 523/5123 and 5123 = 5 is again the centre-of-mass energy squared. For 
this particular simple example we can write the three-particle phase space in D dimensions as 



d(^3 = d(^2d<^ 



unres 5 



d(^u 



r(i-£ 

Integration over the phase space yields 
2 a. _ / 1 



/^ Ne-2 

5123^*^ ■'fS(l —Xi — JC2— .JCs) (.'Cl.'C2-'C3) 



d(^3 



R 



(0) 



—Cf ( ^ + I- + ^ - i^^ 1 S^ / d^2 



n ' \z^ le 4 12 
We see that in the sum the poles cancel and we obtain the finite result 



SI 



(0) 



fd(^2 ( 



2Re^f*^« 



+ / #3 



^ 



(0) 



3^ a, 

-Cf— 
4 71 



d(^2 



SI 



(0) 



(106) 



+ 0(£). (107) 



+ 0(£). (108) 



In this example we have seen the cancellation of the infrared (soft and coUinear) singularities 
between the virtual and the real corrections according to the Kinoshita-Lee-Nauenberg theorem. 
In this example we integrated over the phase space of all final state particles. In practise one is 
often interested in differential distributions. In these cases the cancellation is technically more 
complicated, as the different contributions live on phase spaces of different dimensions and one 
integrates only over restricted regions of phase space. Methods to overcome this obstacle are 
known under the name "phase-space slicing" and "subtraction method" [67-74]. 

The Kinoshita-Lee-Nauenberg theorem is related to the finite experimental resolution in de- 
tecting final state particles. In addition we have to discuss initial state particles. Let us go back 
to eq. (14). The differential cross section we can write schematically 

dOHiH2 = Y^ I dxifHi^a{xi) I dX2fH2-^b{x2)d(5ab{xi,X2), (109) 
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where fn^aix) is the parton distribution function, giving us the probability to find a parton 
of type a in a hadron of type H carrying a fraction x to x + dx of the hadron's momentum. 
d(5abixi,X2) is the differential cross section for the scattering of partons a and b. Now let us look 
at the parton distribution function fa^b of a parton inside another parton. At leading order this 
function is trivially given by 5a^5(l —x), but already at the next order a parton can radiate off 
another parton and thus loose some of its momentum and/or convert to another flavour. One finds 
in D dimensions 

1 a 

fa^b{x,e) = bab8{l-x)--^P^^,{x) + 0{a^.), (110) 

where P^^i, is the lowest order Altarelli-Parisi splitting function. To calculate a cross section 
dOHiHi at NLO involving parton densities one first calculates the cross section d&ab where the 
hadrons Hi and H2 are replaced by partons a and b to NLO: 

d^ab = dol, + ^dolt, + 0{a]) (111) 

The hard scattering part d(5ab is then obtained by inserting the perturbative expansions for d&ab 
and fa^b iiito the factorisation formula. 

dolb + ^dol = do:, + ^dol-l^l^ldxiP'^^,da%-l^l^ldx2P',^,dcyl. 
One therefore obtains for the LO- and the NLO-terms of the hard scattering part 

doib = do\, + lY. I d^iPa-^cdo'l. + lY. I d^iPUdd^li- (112) 

^ c J ^ d ■^ 

The last two terms remove the coUinear initial state singularities in (io^^. 

7 Feynman integrals and periods 

In the previous section we have seen how all divergences disappear in the final result. However 
in intermediate steps of a calculation we will in general have to deal with expressions which 
contain poles in the regularisation parameter £. Let us go back to our general Feynman integral 
as in eq. (81). We multiply this integral with e'l'^^, which avoids the occurrence of Euler's 
constant in the final result: 

Ig = e'^^'^ ^ / d"xd{\-Y,Xi)[[[dXjx/ — -— — ^. (113) 

nr(v,) io -1 Vi=i / ^ 

7=1 
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This integral has a Laurent series in £. For a graph with / loops the highest pole of the corre- 
sponding Laurent series is of power (2/): 

j=-2l 

We see that there are three possibilities how poles in £ can arise from the integral in eq. (113): 

First of all the Gamma-function r(v — lD/2) of the prefactor can give rise to a (single) pole 
if the argument of this function is close to zero or to a negative integer value. This divergence is 
called the overall ultraviolet divergence. 

Secondly, we consider the polynomial Zl . Depending on the exponent V — (/ + 1 )D/2 of 11 the 
vanishing of the polynomial U in some part of the integration region can lead to poles in £ after 
integration. From the definition of U in eq. (91) one sees that each term of the expanded form 
of the polynomial 11 has coefficient + 1 , therefore U can only vanish if some of the Feynman 
parameters are equal to zero. In other words, U is non-zero (and positive) inside the integration 
region, but may vanish on the boundary of the integration region. Poles in £ resulting from the 
vanishing of 11 are related to ultraviolet sub-divergences. 

Thirdly, we consider the polynomial ^ . In an analytic calculation one often considers the 
Feynman integral in the Euclidean region. The Euclidean region is defined as the region, where 
all invariants (pt^ + pi^ + ... + Pii^' are negative or zero, and all internal masses are positive 
or zero. The result in the physical region is then obtained by analytic continuation. It can be 
shown that in the Euclidean region the polynomial y^ is also non-zero (and positive) inside the 
integration region. Therefore under the assumption that the external kinematics is within the 
Euclidean region the polynomial J can only vanish on the boundary of the integration region, 
similar to what has been observed for the the polynomial U . Depending on the exponent v — ID/ 2 
of 5 the vanishing of the polynomial J on the boundary of the integration region may lead to 
poles in £ after integration. These poles are related to infrared divergences. 

Now let us consider the integral in the Euclidean region and let us further assume that all 
values of kinematical invariants and masses are given by rational numbers. Then it can shown 
that all coefficients Cj in eq. (114) are periods [47]. I should first say what a period actually is: 
There are several equivalent definitions for a period, but probably the most accessible definition 
is the following [75]: A period is a complex number whose real and imaginary parts are values 
of absolutely convergent integrals of rational functions with rational coefficients, over domains 
in M" given by polynomial inequalities with rational coefficients. The number of periods is a 
countable set. Any rational and algebraic number is a period, but there are also transcendental 
numbers, which are periods. An example is the number 71, which can be expressed through the 
integral 

n = ff dxdy. (115) 

x2+3;2<l 

The integral on the r.h.s. clearly shows that ;r is a period. On the other hand, it is conjectured that 
the basis of the natural logarithm e and Euler's constant y^ are not periods. Although there are 
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uncountably many numbers, which are not periods, only very recently an example for a number 
which is not a period has been found [76]. 

The proof that all coefficients in eq. (114) are periods is constructive [47] and based on sector 
decomposition [77-82]. The method can be used to compute numerically each coefficient of the 
Laurent expansion. This is a very reliable method, but unfortunately also a little bit slow. 

8 Shuffle algebras 

Before we continue the discussion of loop integrals, it is useful to discuss first shuffle algebras 
and generalisations thereof from an algebraic viewpoint. Consider a set of letters A. The set A is 
called the alphabet. A word is an ordered sequence of letters: 

w = hh-Jk- (116) 

The word of length zero is denoted by e. Let ^ be a field and consider the vector space of words 
over K. A shuffle algebra JA on the vector space of words is defined by 

{hh---k)-ik+l---lr) = Y, ^o(l)^o(2)---^o(r)5 (117) 

shuffles a 

where the sum runs over all permutations a, which preserve the relative order of 1, 2, ..., A: and of 
A: + 1 , . . . , r. The name "shuffle algebra" is related to the analogy of shuffling cards: If a deck of 
cards is split into two parts and then shuffled, the relative order within the two individual parts is 
conserved. A shuffle algebra is also known under the name "mould symmetral" [83]. The empty 
word e is the unit in this algebra: 

e-w = w-e = w. (118) 

A recursive definition of the shuffle product is given by 

(hh-lk) ■ (Ik+l-lr) = h [{h-k) ■ {lk+l...lr)]+lk+i [{hh-lk) ■ (4+2.../.)] . (119) 

It is well known fact that the shuffle algebra is actually a (non-cocommutative) Hopf algebra [84]. 
In this context let us briefly review the definitions of a coalgebra, a bialgebra and a Hopf algebra, 
which are closely related: First note that the unit in an algebra can be viewed as a map from K 
to A and that the multiplication can be viewed as a map from the tensor product A (g) A to A (e.g. 
one takes two elements from A, multiplies them and gets one element out). 

A coalgebra has instead of multiplication and unit the dual structures: a comultiplication A 
and a counit e. The counit is a map from A to K, whereas comultiplication is a map from A to A (g) 
A. Note that comultiplication and counit go in the reverse direction compared to multiplication 
and unit. We will always assume that the comultiplication is coassociative. The general form of 
the coproduct is 

A(a) - £«!^^®«r\ (120) 

i 
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where a\ denotes an element of A appearing in the first slot of A (g) A and a\ correspond- 
ingly denotes an element of A appearing in the second slot. Sweedler's notation [85] consists in 
dropping the dummy index / and the summation symbol: 

A(a) = a(i)®a(2) (121) 

The sum is implicitly understood. This is similar to Einstein's summation convention, except 
that the dummy summation index / is also dropped. The superscripts '^^^ and *^^) indicate that a 
sum is involved. 

A bialgebra is an algebra and a coalgebra at the same time, such that the two structures are 
compatible with each other. Using Sweedler's notation, the compatibility between the multipli- 
cation and comultiplication is expressed as 

^{a■b) = U^^-b^^A<^U'^^-b^^A. {Ill) 

A Hopf algebra is a bialgebra with an additional map from A to A, called the antipode S , 
which fulfils 

a^^^-s{cP-^\ =sU^A-a^'^^ = e-e{a). (123) 

With this background at hand we can now state the coproduct, the counit and the antipode 
for the shuffle algebra: The counit e is given by: 

e(e) = l, e{hl2...ln)=0. (124) 

The coproduct A is given by: 

k 
A{hl2...lk) = £(/,+i...4)®(/i.../,). (125) 

j=o 

The antipode S is given by: 

S{hh...lk) = {-\f hh-i-hh. (126) 

The shuffle algebra is generated by the Lyndon words. If one introduces a lexicographic ordering 
on the letters of the alphabet A, a Lyndon word is defined by the property 

w<v {111) 

for any sub-words u and v such that w = uv. 

An important example for a shuffle algebra are iterated integrals. Let [a,b\ be a segment of 
the real line and fi, fi, ... functions on this interval. Let us define the following iterated integrals: 

b h tk^x 

I{fij2,...Jk,a,b) = J dtifiih) J dt2f2{t2)... J dtkfkih) (128) 

a a a 
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Figure 4: Sketch of the proof for the shuffle product of two iterated integrals. The integral over 
the square is replaced by two integrals over the upper and lower triangle. 



For fixed a and b we have a shuffle algebra: 

I{fhf2,-,fk;a,b)-I{fk+i,...,fr;a,b) 



E ^(/cT(l)./a(2),-,/o(r);«,^), (129) 



shuffles O 



where the sum runs over all permutations o, which preserve the relative order of 1,2, ...,fc and 
of fc+ 1, ...,r. The proof is sketched in fig. 4. The two outermost integrations are recursively 
replaced by integrations over the upper and lower triangle. 

We now consider generalisations of shuffle algebras. Assume that for the set of letters we 
have an additional operation 



A(g)A-^A, 



(130) 



which is commutative and associative. Then we can define a new product of words recursively 
through 



{hh-'-h)* {h+\---lr 



h [{h-h) * {lk+\---lr)]+lk+\ [{hh-h) * [h+i-lr)] 

+ ihJk+l)[il2...lk)Hlk+2...lr)]. (131) 



This product is a generalisation of the shuffle product and differs from the recursive definition of 
the shuffle product in eq. (1 19) through the extra term in the last line. This modified product is 
known under the names quasi-shuffle product [86], mixable shuffle product [87], stuffle product 
[88] or mould Symmetrel [83]. Quasi-shuffle algebras are Hopf algebras. Comultiplication and 
counit are defined as for the shuffle algebras. The counit e is given by: 



The coproduct A is given by: 

A{hl2...lk) = 
The antipode S is recursively defined through 



e{hl2...ln)=0. 



Y^{ij+i...k)^{h...ij) 



k-l 



(132) 



(133) 



s {hh-h) = -/1/2-4- £5 (/;+i---4) * {h-lj) . 

7=1 



(134) 
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Figure 5: Sketch of the proof for the quasi-shuffle product of nested sums. The sum over the 
square is replaced by the sum over the three regions on the r.h.s. 

An example for a quasi-shuffle algebra are nested sums. Let n^ and n^ be integers with n^ < n^ 
and let fi, fi, ••• be functions defined on the integers. We consider the following nested sums: 



S{fi,f2,---,fk',na,nh) 



nt 



/,-l 



I Mh) I /2O2)../ £ fk{ik) 

\=na ii=na ik=na 



(135) 



For fixed Ua and n^ we have a quasi-shuffle algebra: 

S{fij2,---Jk\na,nb)*S{fk+i,...Jr\na,ni,) = 

Y, h(h)S{f2,...Jk,na,ii-l)*S{fk+i,...Jr;naJi-l] 



"b 



+ Y, fk{Jl)S{fij2,...Jk\naJi-l)*S{fk-^2,...Jr\naJl-l] 
rib 

+ L fdi)fk{i) Sif2,---Jk;na,i - I) * S{fk+2, ...Jr;na,i - I) 



(136) 



Note that the product of two letters corresponds to the point- wise product of the two functions: 

(/n/,)W = Mn)fjin). (137) 

The proof that nested sums obey the quasi-shuffle algebra is sketched in Fig. 5. The outermost 
sums of the nested sums on the l.h.s of (136) are split into the three regions indicated in Fig. 5. 



9 Multiple polylogarithms 

In the previous section we have seen that iterated integrals form a shuffle algebra, while nested 
sums form a quasi-shuffle algebra. In this context multiple polylogarithms form an interesting 
class of functions. They have a representation as iterated integrals as well as nested sums. There- 
fore multiple polylogarithms form a shuffle algebra as well as a quasi-shuffle algebra. The two 
algebra structures are independent. Let us start with the representation as nested sums. The 
multiple polylogarithms are defined by [89-92] 



Li, 



mi,...,mi, 



{xu...,Xk) 



ii>i2>--->'k>0 



H' 



ir^ 



(138) 
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The multiple polylogarithms are generalisations of the classical polylogarithms Li„(ji:), whose 
most prominent examples are 



OO 7, OO 71 



Ui{x) = £ — = -ln(l -x), U2{x) = £ — , (139) 



;i = l '1 ii = l h 



as well as Nielsen's generalised polylogarithms [93] 



Sn,p{x) = Li„+i,i,...,i(.x:, 1,..^, (140) 

p-l 

and the harmonic polylogarithms [94, 95] 

^mi,...,mii\X) ^ i^lmi,...,mj.(,-f7 J^5 •••? -U • (J^^J^) 

yt-1 

In addition, multiple polylogarithms have an integral representation. To discuss the integral 
representation it is convenient to introduce for Zky^O the following functions 

GUu...,z,y) = } Jl^'j Jl^.y J!^. (142) 

J h-Zl J t2-Z2 J tk-Zk 


In this definition one variable is redundant due to the following scaling relation: 

G{zu...,Zk,y) = G{xzi,...,xzk;xy) (143) 

If one further defines 

g{z;y) = -^, (144) 

y-z 



then one has 



—G{zu...,zk;y) = g{zi;y)G{z2,...,zk;y) (145) 



and 



G{zuZ2,---,Zk;y) = dt g{zi;t)G{z2,---,Zk;t). (146) 


One can slightly enlarge the set and define G(0, ..., 0; j) with k zeros for zi to Zk to be 

G{o,...Ay) = ^(M'- (147) 
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This permits us to allow trailing zeros in the sequence (zi , ..., Zyt) by defining the function G with 
trailing zeros via (146) and (147). To relate the multiple polylogarithms to the functions G it is 
convenient to introduce the following short-hand notation: 

Gmi,-,mk{zu...,Zk;y) = G(0,...,0,zi,...,Zfe_i,0...,0,z^;y) (148) 

mi — I rnk — l 

Here, all Zj for j = 1, ...,fc are assumed to be non-zero. One then finds 

\Xl X\X2 Xi...Xk J 

The inverse formula reads 



r ( • \ - (-U^^\ ( 1. ^ Zk-i \ 

\Z\ Zl Zk J 



(150) 



Eq. (149) together with (148) and (142) defines an integral representation for the multiple poly- 
logarithms. 

Up to now we treated multiple polylogarithms from an algebraic point of view. Equally 
important are the analytical properties, which are needed for an efficient numerical evaluation. 
As an example I first discuss the numerical evaluation of the dilogarithm [96]: 

Li2(..) = -/j.i^lii^=f;^ (151) 

J t , n^ 

The power series expansion can be evaluated numerically, provided |x| < 1. Using the functional 
equations 

Li2W = -Li2Q)-y-^(ln(-^))^ 

Li2(;c) = -Li2(l-;c) + — -ln(jc)ln(l-;c). (152) 

6 

any argument of the dilogarithm can be mapped into the region \x\ < 1 and — 1 < Kq{x) < 1/2. 
The numerical computation can be accelerated by using an expansion in [— ln(l — x)] and the 
Bernoulli numbers 5,: 

CO n 

Li2W = l77-W(-Ml--^))'+^ (153) 

,=0 I' + ^J- 

The generalisation to multiple polylogarithms proceeds along the same lines [97]: Using the 
integral representation eq. (142) one transforms all arguments into a region, where one has a 
converging power series expansion. In this region eq. (138) may be used. However it is advanta- 
geous to speed up the convergence of the power series expansion. This is done as follows: The 
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multiple polylogarithms satisfy the Holder convolution [88]. For zi 7^ 1 and Zw 7^ this identity 
reads 

G(zi,...,z,;l)= (154) 

52 (-1)^g( 1-z;, l-zj_i,...,l-zi;l j g(zj+i,...,z^;- j . 

The Holder convolution can be used to accelerate the convergence for the series representation 
of the multiple polylogarithms. 

10 From Feynman integrals to multiple polylogarithms 

In sect. 5 we saw that the Feynman parameter integrals depend on two graph polynomials fl 
and jF , which are homogeneous functions of the Feynman parameters. In this section we will 
discuss how multiple polylogarithms arise in the calculation of Feynman parameter integrals. 
We will discuss two approaches. In the first approach one uses a Mellin-Bames transformation 
and sums up residues. This leads to the sum representation of multiple polylogarithms. In the 
second approach one first derives a differential equation for the Feynman parameter integral, 
which is then solved by an ansatz in terms of the iterated integral representation of multiple 
polylogarithms. 

Let us start with the first approach. Assume for the moment that the two graph polynomials 
11 and ^ are absent from the Feynman parameter integral. In this case we have 

1 / \ n r(Vj) 

" fid.,//-'] 6(1 -t^.) ^ ^£^--. (155) 







j=i / i=i r(vi + ...+v„) 



With the help of the Mellin-Barnes transformation we now reduce the general case to eq. (155). 
The Mellin-Bames transformation reads 



100 ;oo 

1 1 



(Ai+A2 + ...+A„) ^ = :^,^— — — j- / dOi... I dOn-i (156) 



r(c) (271/) 
xr(-oi)...r(-o,-i)r(oi + ...+o„-i+c)A;^'...A;;riX 






Each contour is such that the poles of r(— o) are to the right and the poles of r(o + c) are to 
the left. This transformation can be used to convert the sum of monomials of the polynomials 11 
and IT into a product, such that all Feynman parameter integrals are of the form of eq. (155). As 
this transformation converts sums into products it is the "inverse" of Feynman parametrisation. 
Eq. (156) is derived from the theory of Mellin transformations: Let h{x) be a function which is 
bounded by a power law for x — )> and jc — )> 00, e.g. 

\h{x)\<Kx^'^ forjc^O, 
\h(x)\<K'x'' for x^ 00. (157) 
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Then the Mellin transform is defined for cq < Re o < ci by 



h^{o) = dxh{x)x''-K (158) 


The inverse Mellin transform is given by 

y+ioo 

h{x) = — I dahlia) x-"". (159) 

Y-(oo 

The integration contour is parallel to the imaginary axis and cq < Re y < ci. As an example for 
the Mellin transform we consider the function 

h{x) = 7-^^ (160) 

with Mellin transform /z 51^ (a) = r(-a)r(o + c)/r(c). ForRe(-c) <ReY<Owehave 

^ _ 1 7",^r(-a)r(a+.)^_„ ^,^,^ 



{l+xY 2ni J r(c) 

From eq. (161) one obtains with x = B/A the Mellin-Barnes formula 

y+ioa 

(A+Br = — /jaIt^P^A-5— . (162) 

r(c) 

Eq. (156) is then obtained by repeated use of eq. (162). 

With the help of eq. (155) and eq. (156) we may exchange the Feynman parameter integrals 
against multiple contour integrals. A single contour integral is of the form 

y+ioo 



1 

2ni „ 


/ do 


Y- 


loo 



2ni J ^ r(a + C2)...r(o + Cp) T{-<5 + di)..T{-<5 + dq) ^ 

y—ioo 

If max (Re(—(3i ),..., Re(—flm)) < min(Re(Z?i), ...,Re(Z?„)) the contour can be chosen as a straight 
line parallel to the imaginary axis with 

max(Re(-ai),...,Re(-am)) < Rey < min (Re(Z?i ),..., Re(Z?„) ) , (164) 

otherwise the contour is indented, such that the residues of r(o + ai), ..., T{(5 + am) are to the 
right of the contour, whereas the residues of Y{—(5 + b\), ..., T{—(5 + bn) are to the left of the 
contour. The integral eq. (163) is most conveniently evaluated with the help of the residuum 
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theorem by closing the contour to the left or to the right. To sum up all residues which lie inside 
the contour it is useful to know the residues of the Gamma function: 

(-1)" (-1)" 

res {r{o + a),o= -a-n) = - — p-, res (r(-a + a),o = a + n) = -^ — p-. (165) 

In general there are multiple contour integrals, and as a consequence one obtains multiple sums. 
In particular simple cases the contour integrals can be performed in closed form with the help of 
two lemmas of Barnes. Barnes first lemma states that 

if none of the poles of r{a + a)r{b + a) coincides with the ones from r{c — a)r{d — c). Barnes 
second lemma reads 

ioo 

1 f , r(a+o)r(b+o)r(c+o)r(d-o)r(e-o) 

- — / do 



-loo 



Ini J T{a + b + c + d + e + <3) 

o 

T{a + d)T{b + d)T{c + d)T{a + e)T{b + e)T{c + e) 
T{a + b + d + e)T{a + c + d + e)T{b + c + d + e) 



(167) 



Although the Mellin-Bames transformation has been known for a long time, the method has seen 
a revival in applications in recent years [46, 98-1 12]. 

Having collected all residues, one obtains multiple sums. The task is then to expand all 
terms in the dimensional regularisation parameter £ and to re-express the resulting multiple sums 
in terms of known functions. It depends on the form of the multiple sums if this can be done 
systematically. The following types of multiple sums occur often and can be evaluated further 
systematically: 
Type A: 

y r(/ + ai) r(/ + at) i 

rorO'+«i)"TO'+<)'' 

Up to prefactors the hyper-geometric functions y+iFy fall into this class. 
Type B: 



oo oo 



■SJ 



^ ^ r(/ + ai) T{i + ak)T{j + bi) r{j + bi)r{i + j + ci) Tji + j + c^j ,■ 

toh^ii+<y"ni+<)nj+b\y"nj+b'i)r{i+j+c\rT{i+j+c'„^''^ 

An example for a function of this type is given by the first Appell function Fi . 
Type C: 



toh^ j Jni + a[y-ni + a',)r{i + j + c[yT{i + j + c'j''^ 
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Here, an example is given by the Kampe de Feriet function Si. 
Type D: 



An example for a function of this type is the second Appell function F2. 

Note that in these examples there are always as many Gamma functions in the numerator as in 
the denominator. We assume that all a„, a'^, bn, b'^, Cn and c', are of the form "integer + const ■ £". 
The generalisation towards the form "rational number + const ■£" is discussed in [113]. The 
task is now to expand these functions systematically into a Laurent series in £. We start with the 
formula for the expansion of the Gamma-function: 

r(n + £)= (168) 

r(l +£)r(n) [1 +£Zi(n - l)+e^Zn{n - 1) +e^Zni{n - 1) + ... + e"-^Zn..A{n - 1)] , 

where Z^^ . ,„j(n) are Euler-Zagier sums defined by 

Zmi,...,m,{n) = E Jli^----n^- (169) 

n>ii>i2>...>ik>0^ k 

This motivates the following definition of a special form of nested sums, called Z-sums [113- 
116]: 



Z{n\mi,...,mk\xi,...,Xk) = V t^---^1^- (1^0) 



- ^1 ^k 

■mi ■■■ : nik' 
n>ii>i2>...>ik>0^ ^ 



k is called the depth of the Z-sum and w = mi + ... + m^t is called the weight. If the sums go to 
infinity (n = 00) the Z-sums are multiple poly logarithms: 

Z{oo-^mi,...,mk;xi,...,Xk) = U^^,,„^mkixh---,Xk). (171) 

For xi = ... = Xk= I the definition reduces to the Euler-Zagier sums [1 17-121]: 

Z(n;mi,...,m;t;l, •••,!) = Z„j,...,„^(n). (172) 

For n — 00 and xi — ... — Xk= I the sum is a multiple ^-value [88, 122]: 

Z(oo;mi,...,mfc;l,...,l) = ^mi,...,™^. (173) 

The usefulness of the Z-sums lies in the fact, that they interpolate between multiple polylog- 
arithms and Euler-Zagier sums. The Z-sums form a quasi-shuffle algebra. In this approach 
multiple poly logarithms appear through eq. (171). 

Let us look as an example again at eq. (49). Setting D = 4 — 2£ we obtain: 

f d^-^'^h 1 I y^(_ .-i-e r(-£)r(i-£) - r(n+£) i 

J in^-^ i-kj) {-kj) i-q) ^ ''''> r(i-2£) hnn+i)^ ""^ ' 

(174) 
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with X = (— ^i2)/(— ■^123)- The simplest way to arrive at the sum representation is to use the 
following Feynman parametrisation: 



1 1 

/ = {-si23y^'^r{l+e) fda fdbb-^-\l-b)-^l-a{l-x)] 



l-e 

uu u [^i—U) [1— ai^i— Ajj 





(175) 



One then expands [1 — a (1 — jc)] ^ according to 

We continue with eq. (174): Expanding r(n + e) according to eq. (168) one obtains: 
^ r(-e)r(l-e)r(l + e)(-..,3)--£^.-. 

r(i-2£) i-x ^^ e::zi 

n 

In this special case all harmonic polylogarithms can be expressed in terms of powers of the 
standard logarithm: 

^l,...,l(l--^) = ^(M". (177) 



This particular example is very simple and one recovers the well-known all-order result 

ni-e)^r{l+e){-sj,,)-'-'l-x-' 

r(i-2£) i^ r^T' ^^^^^ 

which (for this simple example) can also be obtained by direct integration. If we expand this 
result in £ we recover eq. (78). 

An alternative approach to the computation of Feynman parameter integrals is based on dif- 
ferential equations [95, 123-128]. To evaluate these integrals within this approach one first finds 
for each master integral a differential equation, which this master integral has to satisfy. The 
derivative is taken with respect to an external scale, or a ratio of two scales. An example for a 
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one-loop four-point function is given by 



ds 



123 




2{D-3) 



(^123 - ^12) (^123 -512-523) 



2(D-3) 



(5123 - 523) (5123 - 512 - 523) 




The two-point functions on the r.h.s are simpler and can be considered to be known. This equa- 
tion is solved iteratively by an ansatz for the solution as a Laurent expression in £. Each term in 
this Laurent series is a sum of terms, consisting of basis functions times some unknown (and to be 
determined) coefficients. This ansatz is inserted into the differential equation and the unknown 
coefficients are determined order by order from the differential equation. The basis functions are 
taken as a subset of multiple poly logarithms. In this approach the iterated integral representa- 
tion of multiple polylogarithms is the most convenient form. This is immediately clear from the 
simple formula for the derivative as in eq. (145). 

11 Conclusions 

In these lectures I discussed Feynman integrals. After an introduction into the basic techniques, 
the lectures focused on the computation of Feynman parameter integrals, with an emphasis on 
the mathematical structures underlying these computations. One encounters iterated structures 
as nested sums or iterated integrals, which form a Hopf algebra with a shuffle or quasi-shuffle 
product. Of particular importance are multiple polylogarithms. The algebraic properties of these 
functions are very rich: They form at the same time a shuffle algebra as well as a quasi-shuffle 
algebra. Based on these algebraic structures I discussed algorithms which evaluate Feynman 
integrals to multiple polylogarithms. 
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